QCD correction to single top quark production at the ILC 
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\ Single top quarlc production at the ILC can be used to obtain a high precision 



measurements of the the Vtb CKM matrix element as well as the effective thW cou- 
pling. We have calculated the QCD correction for the cross section in the context of 
p an effective vector boson approximation. Our results show a ~ 10% increase due to 

the strong interaction. 
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I. INTRODUCTION 

The top quark stands out as the heaviest known elementary particle and its properties and 
interactions are one of the most important measurements for present and future high energy 
^ . colliderspj. At the Tevatron and at the LHC the process of single Top quark production 
^ ■ has been extensively studied Q]. 

The top quark is likely to provide us with the first clues of physics beyond the Standard 
Model In fact, new physics effects are probal 
backward asymmetry observed at the Tevatron 

The planned International Linear Collider (ILC) will collide electron and positron beams 
at an initial energy of 500 GeV and higher. It will provide a clean environment for the study 
of precision measurements. 

The single top production processes at lepton and photon (e"'"e~, e~e~, 76 and 77) 
colliders have been extensively studied at tree level in Ref. jsj. The reaction 76" — t- ihve^ is 
particularly suitable for precision studies, as it does not have the ti background. Compared 
to the ILC e+e" — tbe~Ue process the 76" reaction can yield a larger production rate and is 



Dly already manifest in the recent forward- 

la. 
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directly proportional to the Vtb term. Further studies, have thus been done for this reaction. 
In particular, the QCD corrections have been studied in Ref. ^|. Their conclusion is that 
the QCD correction is not very large (~ 5%) so that this mode remains very well suited 
for a precise measurement of Vtb- The approach used by [9|] is to use the effective vector 



boson approximation, also known as effective W- approximation jo] (EWA) and to compute 
the QCD loop corrections for the W^'y — )■ tb fusion process. Then, the convolution with 
the fw+/e+{x) distribution function is applied to obtain the correction to the actual e"*"7 



process. We would like to point out that the authors in Ref. {q] have made a very clear 
and thorough presentation of the calculation. In this work we use their analysis on the 
W'^'-f — 7- tb process to estimate the QCD correction for the e~^e~ — ?■ tbe^Ue process of the 
ILC. Here, in addition to the convolution with the boson we will use the effective photon 
(as well as the effective Z-boson) approximation to obtain the QCD correction. We will use 
the same input values for masses and coupling constants, except for the masses of top and 
bottom quarks we take = 173 GeV and nib = 4.2 GeV. 

II. VECTOR BOSON CONTRIBUTIONS AT TREE LEVEL 

At tree level there are 20 diagrams for the e^e^ — t- tbe^Ve processlS]. We can list them 
in three different types: (a) vector boson fusion, (b) vector boson exchange and (c) e^e~ 
annihilation (see Figur^. For the energy range we consider one of the diagrams actually 
corresponds to ti production, where one of the tops decays leptonically. In order to exclude 
tt production from the single top process we discard all events where the invariant mass of 
the decay products (e~,z/e,6) falls inside an interval around the top mass mt — AM < M^^b ^ 
mt + AM. We take the value AM = 20GeV as in Ref. Is]. 

The effective-W approximation relies on the fact that the vector fusion diagrams become 
dominant when heavy particles are produced at very high energy collisions [si]. In general, 
3 conditions should be met for the EWA to work well: (1) The mass of the vector boson 
{My/ or Mz) should be much smaller than its energy, and this can be met if we require 
My <^ (2) for qq production irtq ^ My, this is true for the top quark but not for the 

bottom quark, and (3) One polarization mode should be dominant so that interference effects 
can be neglected. Fortunately, in our case the mode W'-f tb dominates for longitudinal 
W, and the modes with the Z boson WZ — )■ tb give even lower contributions. 
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As expected, this method works very well for tt production at high ^/s and to a lesser 
degree for single top, which in our case can be seen as tb production. In Ref. 9j the QCD 
correction to the process e"^7 — )■ ibue was calculated by doing first the QCD correction to 
the W'^'y fusion into tb and then by taking the convolution with an effective coming 
from the initial positron (see Figure [2]). We follow the same approach by doing the one loop 
QCD correction to W'^'y — )■ tb as well as W^Z — )■ tb and then convolute with the effective 
distribution functions for W~^, 7 and Z: 

cr(e+e" tbuee^) = (1) 

/ dXwfw+/e+{xw) rfx^/-y/e-(x^) 0-(W+7 -> t6)(s) 

+ y2 dxwfw+/e+{xw) dxzfz/e-{xz) (^{W^Z ^tb){s) 

Wl,t,Zl,t w •'^z 

Where, = 2Mv/\/s, s = xwx^s and the structure functions can be found in [6|. The 
tree level cross section for the single top production at the ILC is shown in Fig. |3l The 
exact Born level calculation for the e^e^ — ?■ tbe^Ue process is obtained with Calchep 1^ 
and is shown by the solid line. We can see that the prediction of the EWA (dot-dashed 
curve) is in very good agreement with the exact result for center of mass energies above 1.5 
TeV. However, for the energy range of the ILC the EWA values can be significantly lower. 
In particular, for a/s = lOOOGeV there is a 15% difference and for a/s = 500GeV the EWA 
result be about one half of the exact value. 

There is one aspect of the calculation that is worth mentioning. Because of the kinematics 
of the W^Z — )■ tb process, we run into a divergent behavior as we integrate over the Mandel- 
stam variable t (or the polar angle of the outgoing quark). At a certain value of t the massive 
Z boson can actually decay into bb and this makes the bottom quark propagator to hit a pole 
at this value. We were able to avoid this singularity by setting fc| = instead of Mz- This is 
completely justifiable in the context of the EWA. Let's understand more the importance of 
the assumption My <^ y/s/2. In the complete process (like e"'"e~ — )■ W~^*Z* — )■ tbe^Ue) the 
virtual Z gets a space-like momenta /c| < and is always far from the on-shell condition. 
In fact, the EWA works better when the initial state vector boson momentum square is set 
equal to zero: fc| = 0, = (see Ref. [3] for a detailed discussion). Nevertheless, when 
dealing with a process like ti production one may set A;| = M| as this introduces only a 
small error of order rriz/^/s. It is customary to set the external massive and Z on-shell 
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FIG. 1: The three type of diagrams for the e+e — )• tbe process. 
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FIG. 2: The vector boson fusion diagrams for the W^^{Z) — )■ th process. 



for convenience. However, for the single top process the fact that Z is heavy enough to decay 
into hh is prompting us to implement the A;| = condition in order to avoid the divergent 
behavior. Notice that a similar situation does not apply to the W'^ boson as it cannot decay 



into th. Therefore in our study we choose to keep the on-shell condition kl^ = for the 
initial state but impose fc| = for the Z boson. For the case of the we have checked 



that indeed by setting fc^ = we don't find a significant change in the result. 

Below, we will describe the QCD corrections to the W'^'y and W~^Z processes, includ- 
ing the Dipole substraction method of infrared divergencies. We have followed closely the 
analysis done for the W'^'j mode done by Kuhn et.al. in Ref. 9|. 



III. QCD CORRECTION TO THE W+j{Z) tb PROCESS. 



The QCD loop correction to the W^'~f{Z) — > tb process is given by 9 Feynman diagrams 
(see Fig.2 of 9|). The renormalization procedure involves only the quark's wave function 
and mass parameter. Specific formulas can be found in j^j. Concerning the renormalization 
scale dependence we have also set as at the scale /i = ^/s for our numerical calculation 



5 



" 15 
l> 









WlT 


Born exact 




\ yy EWA lolal 


WtT - 
WlZt 



500 1000 _ 1500 2000 2500 



FIG. 3: The contributions from Ty^7 and Z fusion to the e^e — )• the process. The sohd 
hne shows the exact calculation. 

(it becomes a/J under the convolution). The extraction of IR singularities is done with 
the substraction method of the dipole formalism ^Jj]. This method consists of adding and 
substracting a so-called dipole term: 



e=0 



(2) 



a^^^(iy+7 -^th)= [ [(rfa^),=o - {da"" ® dVd^poie) e=o\ + / [da'' + da"" ® I 

Jtbg Jtb 

Where da^ comes from the real emission W^'~f{Z) — )■ tbg process and da^ ® dVdipoie is the 
substracting dipole term that matches point-wise the singularities associated to the soft 
and/or collinear gluon. Both terms are calculated in = 4 dimensions. In the second 
integral the same dipole term has been partially integrated in the gluon phase space and 
then added to the virtual correction da^ . This sum is performed in c? = 4 — 2e dimensions 
(consistent with the dimensional regularization). 

The general formula for the dipole term is found in Eq. (5.16) of The specific 

expression in our case is: 



da" ® dV, 



dipole 



2kq ■ kt 



\Mo{kgt,kb)\' + {t^b} 



(3) 



where 



'gt 
Vgt,b 

kb 



1 - zt{l - ygt,b 



^gt,b 



]} 



h ■ h 



ygt,b 



{kt + kg) ■ kb ' 
^{l + agt,by -al^ Jzb, 



2 ^9 ' 

SXtb 



kg ■ kt 

Vgt,b = 

Xtb 



2zb 



A 



^gt 



Xtb(l- ygt,b) 

P — kb, P = kiv + k. 



'7 ! 



6 



and Aio{kgt,kb) is the Born level W^'j — > tb amplitude with one modification: the final 
state momenta h and kf, have been replaced by k„f and kh respectively. 

The otKe. va.aMe. a. de«„ed as in ^: „ ^ ^ .J, . ^ 1 . . - 

Xfe = 1 + 2;;, - zt, Xtb = I - zt - Zb, Xtb = K^,Zt,Zb), Xgt = \{l,{kg + kt)'^/s,Zb), and 
X{x,y,z) = a/x^ + ?/2 + 2;^ — 2xy — 2xz — 2yz. 

For the real emission correction we have prepared a Fortran program that integrates the 
cross section for the W~^'~f — )■ tbg process along with dipole substraction. As it turns out, the 
substraction term defined by the dipole formalism in the first integral of Eq. (|2]) is actually 
a very good approximation to the real emission cross section in an important part of the tbg 
phase space, so that the numerical results we obtained were very small: about two orders of 
magnitude below the values obtained for the virtual correction. 

The expression for the dipole term in the virtual correction is: 

where AidiW'^'y — )■ tb) is the Born level amplitude in c? = 4 — 2e dimensions (the fiux 
term of the tb phase space integration is understood). The dipole function is given by 
lgt,b = CF[2P'^ + I'gf}] ( also lgb,t = IgtAt ^ b}) , where P''' and I^^^l are given by Eqs. (5.34) 
and (5.35) in |ll|: 

^ Inp _ _ 1 _ 1 

^ 2e 6 ' ^ ^^^bJ i ^ nt ^ Hb 

+2U,{-p) - 2Lz2(l - p) - - pI) - - pI)} 

Ilfb =- + 3 + ln/ii + ln(l-p,)-21n[(l-pb)2-z,]--^ (5) 
" ' e I- Pb 



2 

Xtb 



Pbil -2pb) + zt\n- ^* 



1 - Pb. 

where p^ = {xtb - Xtb) / {xtb + Xtb) , Pt = {xtb- Xtb + 2zt)/{xtb + Xtb + 2zt), and pb = Pt{t ^ b]. 
These formulas also appear in [Qj], except that in their Eq. (4.14) I^fl the constant term 
should not be 5 but 3. 

r 

Concerning the calculation of da^ , the details can be found in Ref. |9i]. We actually 
worked out this same computation before doing the case for the Z boson. As expected from 
the results shown in Fig. [3] the contribution from the W^Z fusion is much smaller than 
the one from W^'y. In fact, we only considered the correction for the polarizations 
longitudinal and Z transversal as the other possibilities are negligible. 
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FIG. 4: The QCD correction from W'^j and W~^Z fusion to the e~^e~ — )• the~i>e process. The soUd 
Une shows the exact Born level calculation. The Born plus QCD correction is shown in the dashed 
line. 

Our results are shown in Fig. |H The QCD correction for the single top production in the 
e^e~ collision process is of order 10% of the Born level cross section. It will be interesting 
to compare this result based on the effective W-approximation with a future more robust 
calculation based on the complete e^e~ — > the'Ve process. 
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